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Abstract 

We demonstrate the impact that present lepton flavour and CP violation data - 
neutrino oscillations, baryon asymmetry of the universe, flavour violations in charged 
lepton decays and lepton electric dipole moments - have on supersymmetric seesaw 
theories by analysing the class of models based on a U{1) flavour symmetry. The 
fermion U{1) charges are determined to a large extent by the data, the flavour pat- 
terns being naturally defined through their choice. The selected models generically 
predict fj, ^ e^y within the reach of the projected experiments, which could be sen- 
sitive enough to exclude the whole class of them. By now, the present sensitivity 
to /i — ^> 67 already provides stringent bounds on the contribution of the seesaw 
couplings to the lepton electric dipole moments. 
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1 Introduction 



The simplest way to provide for light neutrino masses is the seesaw mechanism 
pp. However, even in its basic type I version, the flavour pattern of the Dirac 
couplings Y^, and the Majorana mass matrix Mr remains largely undetermined. 
Indeed, the already very accurate data on the atmospheric and solar mass differences 
and mixing angles cannot discriminate between effective neutrino mass matrices, m,^, 
with different overall patterns. The potential measurements of other parameters 
in rriiy cannot resolve the indeterminacy which is an intrinsic feature of the heavy 
neutrino decoupling process - see e.g. Complementary informations on the 

Y^, and Mr flavour patterns are to be found in the quantum effects involving the 
heavy neutrinos. It is well known that this is the case of the lepton asymmetry of the 
universe [7j and, assuming supersymmetry, of charged lepton dipole transitions, like 
the lepton flavour violating (LFV) decays Ci — >• £^7 and the CP violating (CPV) 
lepton electric dipole moments (EDM) di jHlCni- Without supersymmetry the large 
LFV discovered in neutrino oscillations and the CPV needed for leptogenesis do 
not lead to any measurable effect in the charged lepton interactions [TT], in spite 
of the high sensitivity reached, in particular, in /x — 67 [12] and de [ISl- Anyhow, 
protecting the Fermi scale with respect to the large mass scales involved in seesaw 
models requires a suitable mechanism - such as the supersymmetry framework, which 
is adopted in this work - and the corresponding new physics at scales not far from 
the TeV region. This is basically the region that lepton experiments are already 
testing for LFV and CPV sources. 

The predictions for leptogenesis and charged lepton dipole transitions have very 
different dependences on Y^ and with respect to neutrino oscillations, yet the ex- 
traction of the constraints on the seesaw parameters is not straightforward. Firstly, 
there are more parameters than experimental constraints; secondly, their analysis 
relies on additional assumptions, e.g., the neutrino initial abundance and the super- 
symmetric masses. Therefore one is forced to reduce the number of parameters by 
considering a model or a class of models at a time. Of course, the interest of these 
analyses depends on the naturalness of the assumptions and the robustness of the 
results with respect to variations in the parameters. 

In this paper, we contribute to this exploratory work and show the power and the 
complementarity of the combined constraints from all these experiments by studying 
a popular class of models: the supersymmetric Froggatt-Nielsen ones where the 
seesaw parameters are restricted by a U{1)f flavour symmetry ^3]. Actually, these 
models ^3 EI seem to best balance theoretical simplicity and consistency with 
experiments ^. Their basic ingredients are: the charges of the three lepton doublets, 
£j {i = 1,2,3), and those of the three heavy right handed neutrinos, nf, one small 
parameter, e, associated to the breaking of U{l)p and to a charge —1. The abelian 
flavour symmetry fixes the matrices Y^, and Mn up to 0(1) coefficients and thus 

^See e.g. ^Hl for a comparison with less or more structured U{1)f models. 
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the e dependence of the various lepton observables in terms of the charges ^i and 
Hi. This class of models has at least seven attractive features : (i) there is only one 
small parameter e besides the B — L breaking scale vb-l] (ii) the hierarchies in Mr 
depend only on nf, (in) the ratio between the (non- vanishing) Y^, couplings of any 
right handed neutrino only depends on corresponding ii differences, so , (iv) my has 
hierarchical eigenvalues depending only on ^i differences; (v) the models are natural 
and the hierarchies do not rely on the unknown 0(1) coefficients; (vi) the texture 
(the zero entries) of Yy is determined by supersymmetry (analytic superpotential); 
(vii) the hierarchies in the Yy elements are fixed by the U{l)p invariance of the 
superpotential. It has one (well known) slightly unattractive feature: some amount 
of tuning is needed to preserve the hierarchy between the solar and the atmospheric 
mass differences. 

Of course, many aspects of these models have been studied in the literature, 
mostly at a time when the experimental situation of neutrino oscillations was more 
poorly established, leaving more freedom for the parameters ^i and e. One of the 
first studies on leptogenesis considered models of this kind and showed that the 
final lepton asymmetry does not depend on M2 and M3 but only on the lightest 
eigenvalue. Mi, of ^H]- LFV has been considered in some models of this class, 
and it was noticed that with 0(1) Yukawa couplings in Yy the predictions are close 
to the experimental limits for /i — 67 [2n]-|2Sj- The significant improvement in 
our knowledge of neutrino oscillation parameters motivates our reappraisal of the 
abelian flavour symmetry models. 

We show that, in this framework, the mass eigenvalues Mi of Mr are already 
almost fixed by present experiments and basically the whole class of models might 
be excluded by future ones In addition, we also derive bounds on the lepton 
EDM from the limits on /i — >■ 67 (which do not apply to other seesaw models) to 
conclude that the phases in Yy alone are not enough to radiatively induce lepton 
EDM even at the level of the planned experiments 1^ . 

We first consider models where, from the choice of the U{l)p charges, no texture 
zeros are natural in Yy, in the sense that they are not required by supersymmetry 
analyticity. A characteristic of this whole class of models is that, in the seesaw 
formula, niatm gets roughly equivalent contributions from all heavy neutrinos. This 
fact by itself puts an upper bound 0(5 x 10^'*)GeV on all the Mj. Oscillation 
experiments then ask for £2 = ^3 = ^1 — 1 and e = 0{msoi/rnatm) (actually close 
to the Cabibbo angle, often chosen for e in models for the quark mass hierarchies). 
The models then yield the prediction Ues = 0{msoi/fnatm), to be tested in the near 
future. The heavy neutrino charges, rii remain free and so do the Mj. 

The next step is to show that the interplay between the constraints from leptoge- 
nesis, LFV and EDM, allow to gain significant information on the spectrum of right 
handed neutrinos in this class of models. With the matm matrix elements fixed up to 
0(1) coefficients, leptogenesis is possible in these models only if Mi = O(10^^)GeV. 



2 



The rate of LFV depends linearly on M3, while EDM depend on the product M2M3. 
The limit on M3 from the fi —>■ ej search presently is two orders of magnitude above 
the leptogenesis value for Mi and is on the way to be improved. Thus, this class of 
models has the nice feature that the comparison between theoretical parameters and 
the more restrictive experimental data is one-to-one in practice, giving also some 
predictions for quantities to be better measured in the future. As a matter of fact, 
EDM turn out to be so small that the only constraint on M2 comes from its defini- 
tion. Ml < M2 < M3. This is because there is a relation between lepton EDM and 
LFV lepton decays while the latter are related by the U{1)f symmetry. Therefore 
all these quantities are bounded from the strong limits on /i — 67, providing further 
tests of the models. 

In this class of models, the strong upper limits on the heavy neutrino masses 
are related to the fact that their contributions to matm are comparable. However, 
the heaviest neutrino can be naturally decoupled from the and r doublets if one 
chooses ris low enough so that the corresponding matrix elements in Yi, are zero by 
the analyticity of the supersymmetric Yukawa couplings. This alternative class of 
models is also studied here and we show that M3 gets close to 0{Mgut)- However, 
the somewhat surprising result is that the would-be decoupled heavy neutrinos are 
even more coupled in YJ. As a consequence these models are already generically 
excluded by the /i — 67 data. Barring a complete decoupling of any heavy neutrino 
from all light ones is enough to select a very restricted, hence predictive, class of 
models when the experimental constraints are imposed. 

A comment is in order on the supersymmetric flavour problem in the context of 
flavour theories, in particular those based on abelian groups. In general, the breaking 
of the U{1) p symmetry defining the small number e induces supersymmetry breaking 
components for both D and F auxiliary fields along U{l)p charged directions, which 
in turn generically generate soft terms with LFV misalignment and CP V phases [27] . 
The flavour non-diagonal contributions to the soft slepton masses come out of order 
(£j — £i)e^^~^' and {rij — nj)e"j~"' with respect to the diagonal slepton masses ^. 
They are analogous (up to charge difference factors) to the pattern of the radiative 
corrections discussed in this paper, since their e dependence is basically fixed by the 
[/(I) F broken symmetry too. A fine-tuning between the mostly constrained tree-level 
and radiative contributions looks very unnatural. Hence, one can naturally assume 
that each contribution is separetely restricted and should be separetely cared for. In 
particular, a suppression mechanism for both the D and F supersymmetry breaking 
associated to the U{1) p symmetry breaking, would lessen the tree-level contributions 
without any impact on the radiative ones. 

^The origin of these flavour problems is related to the structure of the Kahler metric and 
the Yukawa interactions discussed in Appendix B. When the U{1)f is anomalous, its breaking is 
induced by a Fayet-Iliopoulos term and the associated D-terms give rise to an inversed hierarchy 
for the scalar masses and potentially large LFV and CPV in the mass matrices in the charged 
lepton flavour basis |27j . 
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The paper is organised as follows. The models are defined in section 2. The 
resulting effective neutrino mass matrix is confronted to the oscillation data con- 
straints in section 3 and the general structure of rrii, is obtained. In section 4, the 
leptogenesis requirements are imposed on the selected models to obtain the order of 
magnitude of Mi. The seesaw induced LFV effects in supersymmetric models are 
briefly reviewed in section 5, with the limits from the BR(// —>■ ej) upper bound 
given as functions of the slepton and gaugino masses. In the case of the models 
without holomorphic zeros in Y,y, one gets upper limits on M3 which are quite low 
and will be further lowered by the future /x — >■ 67 measurements. Section 6 starts 
with a recall of the general seesaw induced EDM formulae written so to make their 
connection to the LFV transitions manifest. Some details are given in appendix A. 
Then, they are applied to the models discussed here where the indirect limits on 
de from n ^ ej supersede those from the de experimental searches. In section 7, 
models with naturallly vanishing couplings in Y,y are considered and analysed along 
the same lines as in the previous sections. Some details are given in appendix B. 



2 The model 

The model is defined in the lepton sector by the following set of U{1)f charges: 
£i, Ci and n^, {i = 1,2,3) for the lepton electroweak doublets, charged and neutral 
singlets, respectively; hd for the electroweak symmetry breaking Higgs fields with 
v.e.v.'s Vu, Vd, respectively; 2q for the field(s) breaking lepton number by two units 
in the Majorana masses and associated with the scale vb-l- The U{1)f symmetry 
breaking field charge is normalized to —1 and the small fiavour symmetry breaking 
parameter, e, is provided by the ratio of its v.e.v. and the cut-off of the flavour 
theory. 

We flrst concentrate on the class of models such that all entries in the leptonic 

mass matrices are non-zero in the U{1)f basis, as any charge imbalance in the cor- 
responding terms in the effective superpotential can be compensated by analyticity- 
preserving insertions of e. Hence, the flavour structure of the matrices in the leptonic 
sector is determined, up to 0(1) coefficients, by the following non- negative powers 
of e: 



where all the elements of the matrices g, h and r are 0(1). The effective neutrino 
mass matrix results in the expression 

vb-l 
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After studying these models in detail we shall extend our analysis to models where 
some powers in (Q) are negative and the corresponding entry vanishes in the U{1)f 
basis. 

Models with non-zero entries have several interesting properties for the purposes 
of our analysis. Firstly, integrating out the Froggatt-Nielsen field preserves the 
analyticity in the effective theory, at least at the lowest order in e and the powers 
of e are given by the charge imbalance in the effective operators. This is evident in 
(0), where m,^ is independent of the charges rij ^ISj and gives no information on the 
heavy right handed neutrino spectrum. Secondly, the structure - in terms of powers 
of e - of the unitary transformations diagonalizing the Yukawa couplings and the 
Majorana mass matrix in is fixed by the differences between the U{1) charges of 
different generations. As a consequence, in the basis where Yg and Mji are diagonal, 
only the elements of g, r (now diagonal) and h change with respect to eq. ((H), but 
they generically remain 0(1) (for a detailed proof see e.g'. JHl)- Since the abelian 
flavour symmetry only fixes the powers of e up to arbitrary 0(1) coefficients, it is 
convenient and generically equivalent to work directly in the basis where Yg and Mr 
are diagonal. Hence, in the following we choose this basis while keeping the same 
notations as above for simplicity. In particular (j2)) is now in the physical charged 
lepton basis used in the fits of neutrino oscillations, also convenient to study charged 
LEV and CPV transitions, of course. 

The heavy right handed neutrinos possess a characteristic property in this class 
of models: the branch ratios for their decays into each light lepton flavour are all 
of the same order of magnitude and proportional to e^^^ . This property is partially 
lost when some matrix elements in vanish by analiticity as discussed in section 



3 Neutrino Oscillations 

Let us start with the comparison of the effective neutrino mass matrix in Q and 
the corresponding parameters as determined by the experiments. In our analysis we 
adopt the 3 — 0" ranges of ref. pSj : 



m@ = ^Jml -ml = (.05 ± .01) eV , m© = ^J ml - mj = (.0091 ± .0005) eV 

.71 < tan^zs < 1-45 , tan^ia = .65 ± .12 , tan ^13 < .21 . (3) 

Of course, this comparison is present in several papers, e.g. [IHIIIHIIIZI, but the recent 
improvements in the experimental determination of neutrino oscillation parameters 
- in particular in the solar sector - now allow for a basically unique solution, that 
we turn to discuss. 
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First of all, the very large atmospheric mixing angle requires ^2 = ^3 = ^- It is 
convenient to introduce the parameter p = A23/'m'^, where A23 is the sub-derminant 
of the i,j = 2, 3 sub-matrix of and rris is its largest eigenvalue. The bound on 



requires ii > i leaving the larger elements in the i,j = 2, 3 sub-matrix of eq. (0) to 
be of 0(m@). We can take £1 = £ -|- 1 up to a trivial re-scaling of e. From these two 
requirements, the flavour structure of rrii, in this class of models is fixed : 



vb-l 



This structure has also to accomodate the ratio mQ/mQ = 0.18 ± 0.05 as well as 
the relatively large but not maximal value of tan^^i2. One has p ~ cos"^ 6127712/1713 
and tan ^12 = 0{e/p). Then, for ~ and ^ m2, by using the experimental 
range for the solar angle one obtains 0(e) = sin26'i2/2m@ ~ .08 ± .01. In the 
spirit of the whole approach, we take advantage of the 0(1) arbitrariness in the 
matrices h and r to fix e at the central value of mQ/mQ, e = 0.18, from now on. 
From (jH), it follows that Ue^ = 0(e), not so far from the present upper bound 

[niiiHi. 

Clearly, p = 0(e) requires a relatively mild tuning of the 0(1) coefficients, in- 
herent to this class of models jTHj, which remains a weak point of the present class 
of models ^. To illustrate it, in fig. ^ the p distribution is shown for randomly 
generated 0(1) complex matrix elements in h and r and also for the small subset 
of the generated models that fulfil all the experimental constraints - the latter 
magnified by a factor of 10 in the figure. The experimentally allowed region of 
cos^ di2mQ/m@ is also shown for comparison. One can estimate from the figure the 
probability for the generated models to satisfy the required tuning in p. Of course, 
the relatively high precision of the neutrino oscillation data further reduces the 
experimentally viable models to the small fraction shown in the figure. 

The heavy Majorana neutrinos all contribute the same order of magnitude to 
each element of mu independently of their masses, Mj. However, there is an upper 
bound on their masses if their couplings Y^, are to remain perturbative. Indeed, from 
eq. (jH) one has 



2 2 

* Mi vb-l d@ 



•^Models with single right handed neutrino dominance, where the suppression of p becomes 
natural, were specifically proposed in order to overcome this problem . 
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Figure 1: Distribution of p obtained from a sample of models with charges £i — 1 = 
£2 = £3 = i, e = 0.18 and matrix elements rj, hij randomly generated among the 
complex numbers with absolute values between 0.3—3. The upper (blue) distribution 
refers to the whole sample, while the lower (red) distribution is obtained by keeping 
only the models in the sample satisfying all the present neutrino oscillation bounds 
0. The latter has been magnified by a factor of 10 for graphical reasons. The 
(yellow) shaded region represents the present range of cos^ ^12^0/771®. 




where d@ = 0(1) is the analog of m@ for the dimensionless matrix {h'^r ^h)ij e'^'^+'^Ji. 
Since by assumption £ + h2 + Ui > 0, we obtain the upper bound for each Mj 



^ = e2('+'"+"'W. < 0(1) , M@ = ^^5xl0i^GeV . (7) 

Another characteristic property following from £1 = £+1, relevant for the predic- 
tions of the next sections, is the fact that all the right-handed neutrinos, u"^, couple 
to the first family doublet £1 by a factor of e less than to the other two family dou- 
blets. This means that e, are expected to be less important for leptogenesis (the 
decays into fi, and r, z/^ are dominant) and also that LFV and CPV transition 
amplitudes involving the electron are roughly reduced by corresponding factors of 
e = m0/m@. 

We now turn to discuss how other experimental data can give informations on 
the flavour structure of right handed neutrinos, starting with the constraints from 
the assumption of leptogenesis. 
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4 Leptogenesis 



One of the very first explicit calculations of leptogenesis with the seesaw ^H] niech- 
anism adopted models of the same class discussed here. It found that the lepton 
asymmetry critically depends on Mi and is quite independent of the right handed 
neutrino mass hierarchy. Of course, the neutrino oscillation parameters were less 
constrained at that time. Here, we update and slightly generalize it and we ex- 
plicitely find the favoured range of Mi, which is used in the next sections. 

In the notations of ref. |30j, adapted to the supersymmetric case |S1E21ES1, the 
resulting expression for the baryon asymmetry of the universe is: 



r]B = -10 "^eiKf 



X — 1 



In 



1 + a; 



im[(y.y j).^i] 



J) 11 



M] 



Ml 



X 



X 



(8) 



where is the efficiency factor. In the present class of models we can take advantage 
of the X > > 1 limit because both conditions that could invalidate this approximation 
are unnatural. Indeed, without fine-tuning, M\ and M^ cannot be degenerate at the 
level required for a resonant enhancement , and M2 and M3 cannot be a pseudo- 
Dirac pair [2^1 • In the strong washout regime, where the effective neutrino mass 
mi = {YyY},')y\v\l Mx is larger than the equilibrium neutrino mass m* ^ 0.8 x 
lO-'^eV |33j, the dependence of on the initial abundance is very small and well 
approximated by the power law jSHlESEni: ~ (1-5 ± 0.7) x 10"^ eV /ihi. Thus, 
successful leptogenesis yields a general lower limit on Mi as a function of mi, as 
displayed, e.g'., in ref. jSOlISS!- 

In the models considered here, defining the 0(1) matrix Hij = '^i^hikh^i^e'^^''^ , 
one gets: 



3 Mim@ ^ Im[i72j 3^, Mi . Hn ^, . 

£1 = 3 5 — 7.77 3~ = ' mi = m@— — = ^(m®) . (9) 



57r vt 



i ^ Hnrjd® 8n M@ ride 



Hence the approximation for in the strong washout regime is valid and one gets 
for the baryon asymmetry of the universe: 



2 ± 1 10-^eV Ml ^ Im[/f|i] n . . 

z m@ Mq '—^ ri^^ r j 
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Figure 2: Scatter plot in the plane Mi/M@, eq. (fTUj). and mi/m@, eq. (jHl), for 
the models (described in the previous section) satisfying all the neutrino oscillation 
experimental data. The value of Mi/M@ brings a theoretical uncertainty of a factor 
1/2 — 2. For the points in red the sine of the leptogenesis phase is > 0.1 while for 
those in blue it is < 0.1. 



to be compared to the experimental value rjB = (6.3 ± 0.3) x 10"^" "H^. Therefore 
one derives, within the uncertainty of a factor 1/2 — 2, the value of Mi/M@ required 
by leptogenesis. This is shown in fig. |21for the randomly generated models selected 
by the neutrino oscillation data as discussed in section |21 The results favour Mi = 
0(10^^ GeV). It is well known that leptogenesis at such scale hides a potential 
gravitino problem which is solved if the gravitino is either stable |37| or heavier 
than 0(2 TeV) jHHl- Notice that, since ei in eq. © turns out to be close to the DI 
upper bound jSH], the above prediction for Mi is actually close to its general lower 
limit for fhi = 0(m@) [HUl . 

Notice also that in Hji the electron term is reduced by as compared to the fi 
and T ones, so that the leptogenesis phase should mostly result from the latter. 



5 LFV 

We now turn to the experimental data that can provide informations on the heaviest 
right handed neutrino, namely, the LFV decays of the charged leptons. In a theory 
based on broken supersymmetry, ii ij'y arise from the possible misalignment 
between the slepton and the charged lepton mass matrices. The expansion of the 
transition amplitude in terms of the relative misalignment of the slepton masses - or 
insertion approximation - is justified by the fact that these LFV rates are expected 
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to be relatively small, as attested by the bounds on — > 67. There is a rich literature 
on the experimental constraints on this misalignment that emphasize the importance 
of the supersymmetric flavour problem in the lepton sector jlOlE]. Here, we make 
use of the results in 

In supersymmetric theories the radiative corrections are wave function renor- 
malisations by hermitean matrices, whose flavour and CP structures depend on the 
Yukawa couplings. In the context of the supersymmetric seesaw models, the radia- 
tively induced misalignment was emphasized a long time ago [S] and more recent 
analyses |21 1201 HH 12^1 UHl Hlj were prompted by the determination of the neutrino 
oscillation parameters. Starting from CP and flavour conserving soft terms and tak- 
ing Mfi and real and diagonal at Mpi, the flavour and CP violations in slepton 
masses are introduced in the RG equations only by heavy neutrino z/^ loops through 
their flavour changing couplings ^Jjfc^i/fcj, 7^ j, until they successively decouple at 
Mfc yielding corresponding factors ln(Mpi/Mfc). Since couples z/^ to the lepton 
doublets, the seesaw radiative corrections mainly affect the doublet slepton masses 
m|jj-. Assuming universality in soft masses at Mpi, the misalignment in m\ at the 
lowest order in the couplings Yy (deflned at Mpi) is ^j: 



where mg and are the universal soft breaking bilinear and trilinear terms at Mp;, 
respectively, and mz, is the average doublet slepton mass. We have isolated the 
contribution from each state, C^p for later convenience. 

In the basis where Y^ and Mp are diagonal, the mass insertion approximation 
reads 



where Fgusy = 0(1) is a function of supersymmetric masses which includes both the 
chargino and neutralino exchange (see e.g. and references therein). With eq. 
flll|) . the resulting general limits on slepton misalignment jl2] translate into limits 
on the \Cij \ as functions of the supersymmetric masses and tan /3 |23| l^ In mSugra 
one can express the dependence on the susy spectrum in terms of gq, tan/? and two 

"'The figures in |23| only take into account the chargino contribution, not dominant in some 
regions of the supersymmetric masses |42) . The neutralino contributions are included in the plots 
ofref. |4j. 





BR{ii ijj) = 10"^ BR(£, ^ ijUjUi) -^tan^ P\5^/'\'^F, 



susy 



(12) 
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masses: either the bino mass Mi and the average singlet charged slepton mass ffiR, 
or the universal scalar and gaugino mass parameters, respectively mo and Mi/2. The 
upper limits on \C^e\ from the present upper bound BR(/i 67) < 10~^^ are shown 
in fig. 121 for tan/3 = 10. The two choices for are meant to give an idea of the 
theoretical uncertainty: a) Oq = (dotted fines) and b) gq = mo + Mi/2 (solid lines). 
For different cfioices of tan/3 and/or better sensitivities tfian tfie present one, it is 
enougfi to rescale according to (fT^ . Tfie same plot can be used to estimate tfie upper 
bounds on Crfj, from r — /i7 by taking into account tfie factor BR(r /iz/z/) ^ 17%. 

In spite of tfie mucfi poorer determination of e at tfie time, seesaw models of 
tfie present class were among tfie first analysed to get predictions for /i — 67 and 
r fij. It was pointed out tfiat for models wfiose cfiarges allow for 0(1) couplings 
in Yi,, tfie predictions for fj, ^ ej are at tfie level of tfie present limits and tfie 
planned improvement [25 would test tfiem [201 1^ (as can also be inferred from 
fig. 121) • Tfie relevance of r — > /i7, r — > 67 was empfiasized as well. Of course, sucfi 
considerations surpass tfiese particular models. A general discussion of LFV decays 
for tfie tfiree possible classes of seesaw models, defined in terms of tfie dominance 
mecfianism in tfie neutrino masses, was presented in ref. [22] ■ 

Actually, wfiat cfiaracterizes tfie present class of models is tfiat all //"^'s contribu- 
tions to nil, are on an equal- footing, as expressed by eq. (0), so tfiat 



'-'ij ~ "■ki"'kj t t m ■ 



Ml 



k 



KM,s...s,.MJl,^^^ Vfc. (13) 



Tfierefore, from tfie experimental limits on \Cij \ one gets directly bounds on tfie fieav- 
iest Majorana mass M3 - unless, for M3 ~ M2, Cf^ and Cf^ strongly interfere, wfiicfi 
amounts to some tuning of tfie free parameters. On tfie otfier fiand, tfie leptogenesis 
constraint on tfie ligfitest mass Mi, provides a lower limit, Cij > e'^""*"''^^ (10~^ — 10~^), 
fience a lower bound on BR(£j — ^ ^jl)- 

We first discuss in some detail tfie bound on M3 wfiicfi follows from ^ — > 67. 
Tfie rigfit panel of fig. |21sfiows tfie scatter plot of |C^g| - eq. flT^ - as a function of 
Ms/Mq for tfie models fulfilling all neutrino oscillation and leptogenesis constraints 
(of course similar plots apply also to k = 1, 2). Since, witfi tfie value of e fixed as in 
section 121 \C^e\ — 0{M^/ M@), tfie experimental upper bounds on \C^e\ and M^/M@ 
are of tfie same order of magnitude: 



Tr<o(ic,,r-^-)- (14) 

Tfiese upper bounds are displayed in tfie left plot of fig. El Witfi tfie present 
experimental sensitivity, tfiey already provide meaningful limits on M3, in spite of 
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Figure 3: The left plot shows the upper bounds on |C^e| for BR(/i — > 67) < 10^^^ 
and tan/? = 10. The plot, adapted from that in includes both the chargino 
and neutralino contributions and assumes mSugra with oq = (dotted line) and 
Oo = rriQ + Ml/2 (solid one). The right plot shows the scatter plot of |C^e|, eq. 
()13|) . as a function of M^/M@ for the models fulfilling all neutrino oscillation and 
leptogenesis constraints, with Mi < 10^^ GeV. We take M3 to span its allowed range, 
namely Mi < M3 < d@r^. 



the uncertainties in the supersymmetric spectrum. For instance, with tan/5 = 10, 
Ms/Mq < 0.3(1) for the mass region below the orange (yellow, resp.) line, some two 
orders of magnitude below the gauge coupling unification scale! With the planned 
three orders of magnitude improvement in BR(yU — > 67), it would become M^/M@ < 
10~^ in that mass regions, which is close to the Mi range selected by leptogenesis, 
where C^e ^ C'^e ~ (10~^ — 10~^). Therefore, for this class of models, there are fair 
perspectives to detect a positive signal in /z — > 67 searches. Alternatively, barring 
destructive interference among the C^g, the future experiments would rule out this 
whole class of models in the region of the parameter space where IC^d"'^' ^ 10~^. 
For tan/5 > 20, it actually coincides with the area below the red lines which also 
corresponds to a sizeable supersymmetric contribution to {g — 2)^ (see and 
references therein). This clearly shows the power of the combined constraints from 
/i — *• 67 and leptogenesis in this class of models. 

As already noticed, the couplings to e are reduced with respect to /i and r by 
a factor 0(e) = 0{mQ/m@), so that C^^ = 0{C^^ = 0{e)C^^. This simple link to 
the rrii, pattern is a characteristic of the model. In particular, the model predicts: 
BR(/i 67) : BR(r ^ /i7 ) : BR(r 67) = : .17 : .ITe^, while the present 
experimental upper bounds are in the ratios 10~^ : 1 : 1, respectively. Hence, in this 
model, from the experimental limit on yU — > 67, the other two LFV transitions are 
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predicted many orders of magnitude below their present experimental limits. 



6 Lepton EDM 

There is a rich literature about the constraints that lepton EDM searches put on the 
CPV pattern of the MSSM and its extensions I46j . in particular, on the slepton 
mass matrices. It is useful to separate them into two kinds: i) lepton flavour con- 
serving (FC) parameters, namely, the phase of the /z-term and phases in the flavour 
diagonal chirality flipping slepton masses, or A-terms; ii) phases in the lepton flavour 
violating (FV) elements of the slepton mass matrices. Very restrictive conditions on 
some of those phases arise from the present limits on de which are to be considerably 
strengthened in the near future [23]. The present limit on ci^ is not very signiflcant 
but a new experiment [2^] could bring relevant results. The flrst conclusion is that 
supersymmetric models must possess a mechanism to effectively inhibit CP phases 
in the soft terms at Mpi. The second one is that the radiative corrections must 
keep the induced CP phases in the soft terms below the experimental limits. In this 
section we concentrate on the constraints on the radiative corrections arising from 
the Yukawa couplings of the supersymmetric seesaw models. 

The main contributions to diS of either kind depend on the slepton mass param- 
eters as follows: 



PC _ e^MjB , . FV _ e^Mi/xtan/?/^ ^^(sRRrr, 

- 327r2|/.Pcos2^^^^'^^^^^) - 327r2|;,pcos2^H. ^ ' 

(15) 

in the basis where Yi and Mr are real and diagonal and with the FC A-terms 
written as Ai = aiYi-. The functions Ib and depend on the slepton, chargino and 
neutralino masses as deflned in ref. The flrst contributions to df^ are second 
order in the insertion approximation. Only the FV elements of the (hermitian) 
matrices 5^^, 6^^, those responsible for LFV decays, contribute in (jl5j) . Hence d^^ 
is linked to r — >• 67 and /i — * 67, d^^ to r ^ yU7 and fi — > 67, rff ^ to r ^ /i7 and 
r ^ 67 (see, e.g., jl21 for a general discussion). Here, these links are analysed anew 
in the framework of the supersymmetric seesaw mechanism as one source of LFV 
and CPV. 

We start with real fi and no CP or lepton flavour violations in the A-terms at 
Mpi. In particular, Ima^ = 0. Then, CP violations are introduced through the RG 
running only by the heavy neutrino couplings Y,y. The decoupling of right handed 
neutrinos yields an effective MSSM, with neutrino masses and CP phases in the 
slepton masses. As discussed in section ISJ the hermitian matrices appear at one 
loop and the \Cfj\ are bounded by the £i ijj decays. However, since C^^ is real. 
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there is no CPV phase to induce di. This means that the radiatively induced CP 
violations start at the two- loop level and, in the leading log development, come from 
the phases of LFV elements in the products Im(^j. C^jCj^), k ^ k' . By obtaining 
the dependence of df^ and df^ on the C^^ one makes more manifest how they 
are constrained by the LFV experiments. Let us first discuss these relations more 
generally, before turning back to the class of models analysed here. 



6.1 General link between EDM and LFV 



Recently, the dominant two-loop contributions to Ima^ and 5 , 6 in the expres- 
sions of the supersymmetric seesaw radiatively induced EDMs have been calculated 
jni Cni EI- Here, the results are rearranged to make clear the general connection 
between LFV and EDM, by expressing di in terms of the C^y The original results 
in jTUj and some detailed expressions are presented in App. A. 

Assuming mSugra boundary conditions at Mpi for simplicity, the result EI 
for df^ is obtained by replacing in (fT3|) : 



J ( \ ^ ln(Mfc/Mfc/) is:^^kn*k'^ Ma\ 

k>k' j^i 

Notice that since Im(^^ Oj) = the FC contribution leads to the sum rule: 



^dJmi = Q. (17) 

i 

To obtain df^, one has to replace in (fT3|) the leading two-loop result (TUlllZ] for 
the double insertion term, which is: 



RR^ xLL\ 



Im((5''^m^(5 



8 {6ml + 2al) {Qml + Sag) 



(An) 



tFV 



— 1—1 

2 

mi 



mt 



k>k 



mi- tan^ (3 I, 



FV 



where mji and m^ are average masses of the sleptons of each chirality, and the 

~fc 

logarithmic dependence on the masses Mk represented by In^, is defined in App. A. 
The sum rule ()17p is not valid as the FV contribution rather leads to the sum rule: 

5^mic/i = 0. (19) 
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There are other three points to notice with respect to these expressions: 

i) The EDM phase resuhs from the (three different) relative phases between C^j and 
C^j, namely between the and i/^, contributions to ii — > ij'y. 

ii) The phase in the quark CKM matrix is large and there is no known reason to 
think of smaller phases in the lepton sector, thus one expects the EDM phases to 
be generically large. 

in) As already discussed, \C^j\ are constrained by the experimental limits on the 
LEV decays ii — *• ij 7. However, at present the only significant bound from LEV 
decays is that on \C^e\ from — > 67. 

It has been shown in ref. ^ that with (UHl) and (HHD, df^ > df^ for tan /3 > 10. 
Indeed, the EC contribution vanishes for ao = 0, while for ao ~ mo, as df^ /df^ oc 
(tan/5)^, the EV contribution dominates over the EC one for large tan/5 in spite of 
their different sparticle mass dependence. We refer to [TOj for general analytical and 
numerical studies of these seesaw mechanism as a source of lepton EDM. 

Eor those reasons, we mainly concentrate on df^. In fig. |3J the limit on /f^ 
obtained from the bound de ^ 10~^^ e cm is shown for tan/5 = 10 and ao = 
(dotted lines) or ao = rriQ + M1/2 (solid ones). With the present bound the limits 
are 10^ larger and quite useless. Indeed, if the Yukawa couplings are at most 0(1) 
and the masses in the range selected by neutrino masses and leptogenesis, one 
has < 0(10^'^). Then, only for light sparticles and large tan/3 the present 
experimental bound on d^ constrains I^^ to be below this upper bound [TU] . 

Einally, from the explicit expression for (fTBj) , 



one sees that the limit on de from LEV depends more on r — >■ 67 whose experimental 
upper bounds are still very far from providing useful constraints. Therefore this 
comparison is only possible in specific models and the class of models studied in this 
paper is particularly predictive in this respect as already discussed in section El 

6.2 EDM in this class of models 

In the class of models discussed in the previous sections, the C^j are all related, 
so that and C^^ are generically 0(e)C^^. Therefore they are all constrained 
by the most significant limit, that from BR(/x — *■ ej) on C^^. The potentially 
most important contribution to dg comes from the first term of I^^ in ()20p and is 
constrained by the decay r 67. Replacing this term with eq. (fT^ one obtains. 




(20) 
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Lines of constant |/f 




100 200 300 400 500 100 200 300 400 500 
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Figure 4: Upper limit on from < 10~^^ e cm with tan (3 = 10, = mo+Mi/2 
(solid line) and Oq = (dotted line). For larger values of tan/5, multiply the upper 
bound for (10/ tan/3)^. 



where the main contribution comes from A; = 3, k' = 2, the other ones being 
reduced by factors 0(Mi/M2) and 0(Mi/M3), respectively. The /i ^ 67 decay 
giving already a limit on M3, in principle, the experimental limits on the electron 
EDM and (j?T|) yields a bound on M2. However, as we now turn to discuss, the 
seesaw contribution to the lepton EDM comes out too small to give any relevant 
limit in the present class of models. 

The maximum values of (PT|) correspond to M2 = OlM^) and are roughly given 

by 

which is maximum near the maximum value of M3. As shown in section for 
tan /3 > 10 and sparticle masses below the TeV, the present limit on — > 67 
requires M3 < M@, which implies I^^ < 10^^. The planned improvement in fi ^ ej 
searches, at the level of BR(yU — > 67) < 10"^'', would imply M^/M@ < 0.03 in the 
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Figure 5: /f^ as a function of M2IM© for M^/M® = 1 and 0.03 respectively, 
with randomly generated models that fulfill the neutrino oscillation and leptogenesis 
experimental constraints with Mi < 10^^ GeV. 



same region of the MSSM parameter space and, correspondingly, I^^ < lO^^^, thus 
de ^ 10~^^ e cm. Comparing these numbers with the plots in fig. 01 one concludes 
that, even with a bound of 10~^^ e cm on de, this radiative contribution is too low 
to be detected. 

The plots in fig. El show the dependence on M2 for the models that satisfy all 
the experimental constraints of the previous sections. The vertical spreading of 
the points corresponds to the random overall phase. It is worth noting that the 
phases in d^ and in ei are uncorrelated. The former is generically proportional to 
Im(/i3^/t33/i2i^23)' while the latter rather depends on a combination of Im(/i22/'-i2 + 
h2zh\^) and hn.(h^2h\2 + ^33^13); as previously discussed in section HI 

Notice also the predictions 



d,, = O ( e-^^ ] de = 0(6000)4 , dr = oi e'^^^ ) 4 = 0(400)4 , (23) 



niem^ 



very different from the naive proportionality between lepton EDM and mass. Of 
course d^ and dr are also predicted to be extremely small. 
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7 Decoupling the heavy neutrinos from m@ 



The fact that heavy Majorana masses end up pushed down by the experimental 
hmits on LFV is due to the main characteristic of the models studied until now: as 
far as the heavy neutrinos commensurately contribute to the value of m@, the only 
way to reduce the LFV from the Yukawa couplings is to decrease their masses. The 
alternative is to decouple one (or two) heavy neutrino from m@. This is naturally 
done in the abelian flavour symmetry approach by the introduction of the so-called 
analytic zeros, imposed by the flavour symmetry and supersymmetry on the analytic 
superpotential. In the broken flavour symmetry phase, the zeros in the effective 
Yukawa interactions are filled by the diagonalisation of the Kahler metric which 
becomes e dependent. Apparently, their values are diminished by powers of e, and 
so might be their contribution to LFV decays. Actually it turns out that, in the 
basis of the physical charged leptons, the couplings of the would-be decoupled state 
must be so much larger than those in the coupled case that the (partial) decoupling 
option is already experimentally excluded. 

In this section, we try and make the expressions more transparent by omitting the 
0(1) matrix elements defined and displayed in the previous ones. Correspondingly, 
we discuss only orders of magnitude and show the results in this form. More detailed 
calculations and definitions are presented in Appendix B, including a brief account 
of the transformations to the basis chosen here with diagonal Mji and Yg. It is also 
shown that the hierarchical structure, i.e., the e dependence of rrii, as defined by 
is stable under the necessary field redefinitions. From the discussion in section El we 
can already already take £2 = ^3 = ^ = — 1, as in eq. ©, as well the relations in 
© and CZD below. 

In order to decouple from m@, we assume the charge ris to be large enough so 
that n^+i+hu < 0. Since this charge imbalance cannot be compensated analytically, 
it implies that Y^^i = {i = 2,3) in the basis where the abelian flavour charge is 
diagonal. If one also chooses + ii + < 0, z/| decouples from all light neutrinos, 
in this basis, and the model becomes of the so-called 3x2 type, which we skip here 
since it has received a lot of attention in the literature j^. As just explained, we 
are led to choose ii = i + 1 from (0). Hence + ii + hu = and Yi^^i = 0(1). 
Then, from ((7j), we get 

M3 = 0(e-^M@) = O(10^^)GeV , (24) 

close to the GUT scale. Since Mk/M^ = 0(6^"*="^"^), one has in the basis of the 
flavour charge eigenstates: 



F,3i = 0(1) , y,3, =0 (j = 2, 3) , Y,kj = 0{e 



5,1 



(fc = l,2). (25) 



In going to the physical basis where M/j (Yg) is diagonal, there are two steps: 
first to put the Kahler metric in canonical form, then to diagonalise (Y^, resp.). 
While in section |21 it has been remarked that this does not change the structure 
of the matrices in the models previously studied, here they fill the analytic zeros 
Yuzj = (j = 2, 3). In this basis these couplings become (see Appendix B) 

where the 0(1) coefficients come from a small rotation and a small Kahler vielbein, 
the first (second) terms are associated to the transformations of the right-handed 
neutrinos (left-handed charged leptons, resp.). Notice that these couplings are large, 
of 0(e), so large that the model is ruled out by LFV decays. Indeed, keeping only 
the leading z/g and z/| contributions, one gets: 

= . O (i^) 1„ ^ + e f 0(1) + O (^) ) 1„ ^ = 0(1) ^ (27) 

This value is above its MSSM upper hmit IC^el^'^S discussed in sectional even for 
quite heavy gauginos. Mi = O(TeV). The only remaining possibility is to assume 
M2 = 0(M@) and an (unnatural) cancellation between the various 0(1) terms that 
result. In this sense models with decoupling are (i) fixed up to 0(1) coefficients and 
(a) basically already excluded by experiments. Remember that planned experiments 
on — > 67 will reduce lO^er'*^' by more than an order of magnitude. 

Finally, let us just comment on the remaining possibility of decoupling both 
and z/g from m@, namely from both V2 and z/3, so that one neutrino remains exactly 
massless while dominates in m@. This could look interesting because, the tuning 
of p examined in section 01 might become natural [5^. However, it is easy to check 
that this configuration leads to mQ/rriQ = O(e^) and 6^12 = 0(e) requiring unnatural 
coefficients to fit the oscillation data. Of course the same constraints as in the case 
of only one decoupling from m@ also apply here and require a tuning between the 
contributions from z/| and z/g to C*e^, which are both of 0(1). Altogether, this last 
case is even less attractive than the previous one. 

8 Conclusions 

In this paper we have emphasized the complementarity of some of the present and 
future experimental data for the determination of the Dirac and Majorana mass 
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matrices of the super symmetric seesaw models. For this purpose, we have chosen the 
simplest abelian flavour models, in which the structure of the Dirac and Majorana 
masses arc fixed, up to 0(1) parameters with arbitrary phases, by only 8 parameters: 
the 6 lepton flavour charges, the small Froggatt-Nielsen parameter e and the B — L 
breaking scale. We have shown that the available data (basically the atmospheric 
and solar oscillation masses and angles, the baryon asymmetry in the universe and 
the // — > 67 decay) are enough to fix the seesaw mass patterns to a large extent. 
These models predict the // — > 67 decay to be within the reach of the planned 
searches. They also predict a very low contribution to d^. which is bounded in this 
class of models by the limits on — > 67 alone. In other words, the lepton flavour 
abelian flavour models are consistent with experiments but only a few choices of the 
parameters are possible and basically all the selected models will be tested by the 
future experiments. 

The models possessing all these properties are theoretically well-defined effective 
theories, the flavour symmetry is the simplest abelian one. The fact that a few 
experimental data of very different nature can complement each other to constrain 
so much the models is certainly non-trivial. Type 1 leptogenesis always requires one 
relatively light right-handed neutrino, but here it has to be at the upper side of 
the allowed region. It is well-known that the n ^ experiments are testing the 
LEV effects in the slepton masses to the level of the radiative corrections, but here 
they immediately translate into a limit on a right handed neutrino mass eigenvalue 
and the whole series of LFV decays is predicted. Namely, the complementarity of 
the experimental constraints is disentangled into simple relations for the physical 
parameters. 

Of course these models are special in many respects. They predict a hierarchical 
spectrum for the light neutrinos with C/gs close to its present bound, choosing one 
of the possible neutrino mass matrix patterns. They explains the large oscillation 
mixing angles by large mixings in the Dirac masses, favouring large LFV. The ra- 
tio between the solar and the atmospheric mass differences is directly related to a 
fundamental parameter, the scale of the flavour symmetry breaking, at the price of 
some tuning. Whenever all the right handed neutrinos couple to all left handed ones 
in the Yukawa matrix, they all contribute to the atmospheric neutrino oscillations. 
If one right state (not a mass eigenstate) decouples from the /i and the r doublets, 
one Majorana mass eigenvalue can be quite large, close to the GUT scale. With 
the similar decoupling of another state, one gets the pattern corresponding to the 
dominance of the atmospheric oscillations by the lightest right handed neutrino, but 
777.0/777.® and 9i2 come out too small. 

Thanks to supersymmetry, in the class of models examined here all these situa- 
tions are naturally realized in the sense that they are ensured by the choice of the 
arbitrary charges. The results are quite independent of the various 0(1) free pa- 
rameters, meaning that the conclusions are technically robust. The dependence on 
the sparticle masses is taken into account by analysing the results for a large range 
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of values of the most important parameters (slepton and gaugino masses) . On this 
basis, our main conclusions here are: (i) models involving decoupling through su- 
persymmetric zeros are basically excluded and (ii) the whole class of models will be 
generically excluded in the near future, by the experimental bounds on the /i — > 67 
decay - once the constraints from the present data on neutrino oscillations and 
leptogenesis are imposed. 

General predictions of these models are the relations among the LFV charged 
lepton decay rates, and the very low upper bounds on the contributions to de- 
Actually the latter are generically expected HHI if the seesaw couplings are the 
only source of CPV. Instead, if the couplings of the GUT heavy colour triplets 
related to proton decay are also introduced, one can achieve much larger values of 
de lini Uni- In SU{5) GUT, with the 5-plet and singlet charges fixed by the pure 
seesaw contributions, the 10-plet ones by the quark masses and mixings, the limits 
on de are then predicted [^0] up to uncertainties in 0(1) coefficients in the seesaw 
model, their phases, and the amounts of non-minimality in the GUT required to fit 
the data (fermion masses and mixings, proton decay). This analysis of de in a GUT 
framework is also supported by the 0{mQ/m@) value of e obtained in the present 
paper, quite close indeed to the Cabibbo angle which is the natural choice in the 
quark sector. 
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Appendix A 

The FC contribution of Ai = aiYi^ has been calculated in refs. jH] and ^U] at lowest 
order^. In order to treat the various thresholds in a compact way, it is convenient 
to introduce the function defined as In^ = \n{Ma/Mb) where a,b = 1,2,3,4 with 
the identification M4 = Mpi, as well as the matrices Na = Y^PaY^ with P3 = I, 
P2 =diag(l, 1, 0), Pi =diag(l, 0, 0), and it is understood that the Yukawas are taken 
at Mpi. Then [fn] : 

Im(a,) = -^ao/f ^ , 1^ = Y,K''K^'^HNaN,h (28) 

■'The coefficients in the papers of |9j numerically disagree with each other and with the result 
in ref. The differences are irrelevant for the qualitative conclusions herein. Ref. 07] find, on 
the contrary, a smaller FC contribution. 
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where the sum is over a, 6 = 1, 2, 3. 



The contribution from two insertions of FV slepton mass matrix elements ob- 
tained in [ini EZj reads: 



tLL\ 



m: 



(An) 



nig — ^ tan P 



{6ml + 2al){6ml + 3al 



— 2 — 2 



OJjFV 



tFV 



^-^ mi 



a>b 



(29) 



where m^j and mi are average R and L slepton masses at low energy, and Yi ~ 
me tan P/mt has been substituted. 

The FC contribution vanishes for = 0. On the other hand, for ~ m^, it 
turns out that the ratio rff^/rff" is well approximated by the simple rule jTUI 



df^ Inf Aan/3\ , , 



df^ ■ 10 V 10 

where fi has been fixed by electroweak symmetry breaking. Thus, the FV contribu- 
tion is in general favoured by a factor (tan/3/10)'^. A detailed comparison of these 
two contributions is in [lOj, together with limits on J^^^'^'-^ from the present and 
planned sensitivity to d^ and d^. 

The above expressions for I^^ and I^^ can be rewritten in terms of the C's 
related to LFV as follows: 



jFV ^ J2 Jf , = K lmiC^^C% , (31) 

z — ^ mi 

a>b ^ 

~n , o Ini^ In? , , / Inf ' In^ \ , , 

K^.„?-2^. (i- ^ (32) 



Hence, explicitely: 



FV 



jFV 
-'2 



jFV 
-'3 



a>b 
a>b 



a>b 



m. 



mt 



mz 



mt 



m, 



--flm(Q3C: 



mt 



mt 



(33) 
(34) 
(35) 
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On the other hand, for the FC contribution, 

= E il"'^^'' ' ^l^'^^"" = ^Mcch ■ (36) 

Exphcitely: 




It is then manifest that is hnked to r ^ 67 and n 67, ci^ to r ^ /i7 and 
/i 67, and dr to r ^ /X7 and r — ej. 

Appendix B 

This appendix contains some details of the results discussed in section ((7)), namely, 
when a heavy neutrino partially decouples from the /i and the r. In the basis where 
U{1)f is diagonal this is naturally implemented by supersymmetry zeros in the Yi^ 
matrix. In the more convenient basis where and M/j are diagonal, these zeros are 
filled by: i) the vielbein redefinitions of both the light and heavy states [HI], and ii) 
the unitary transformations - that remain free in the vielbein - to go into this basis, 
in this order. In the framework of this paper, where generic 0(1) coefficients are free 
parameters and the calculations are carried on at lowest order, these tranformations 
commute. Therefore, one has to consider the effect of the tranformations: 



Mr V^,^IMr^rVr = M , 

n ^ VieRYuiLVL = Y^ . (40) 

where ^l^r satisfy {iS^Y^ = K'^^ , K^j is the Kahler metric, Vr^l are unitary, M and 
Yi are diagonal. The metric is hermitean and invariant under U{1)f, so that the 
transformation pn|) can be generally written at the lowest order as follows. 
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(eflVff,),, = 5,, + r/ffiiel"'-"^l , {^LVL)^J = k, + r^L^.e'^^-'^l , (41) 

where r^i^^x, are 0(1) matrices whose hermitean parts bring the Kahler metric to the 
canonic form and whose anti-hermitean ones diagonahse M/j and Y^, respectively. 
Of course, the non-diagonal terms are not restricted by analyticity and there are no 
vanishing matrix elements, in general. When applied to Yy this transformation fills 
the possible analytic zeros already at the first order, as follows: 

k k 

If rii and are such that Y^i^ is a supersymmetric zero, it becomes (n^, decreasing 
with i,j): 

k<i k<j 

Notice that, for hierarchical M/j, the rotation angles of Vr are: 

(44) 



sjTr^kk Mi - Ml 

and the vielbein elements are of the same orders of magnitude. 

A few comments are in order here. When the heavy neutrinos are integrated out 
to yield the effective m^, the integration preserves analyticity as it does not involve 
the Kahler metric and the results could dependent only on the transformation ^lVl 
in (j4Up . In the basis of the physical charged leptons one has at the lowest orders: 



^ m^^j + VLki^^^'~^''^mykj + m^ikVLkje^^''^"^ (45) 

with nil, given by eq. (0). It is easily checked that rh^ has the same hierarchies 
in terms of e as rrii,, with the 0(1) coefficients reshuffied by ^lVl- Instead, in 
the calculation of the effects of the RG evolution in LFV and CPV, the metric 
dependence is explicit in the results, even though they remain independent of Vr. 
For leptogenesis, the lepton asymmetry clearly depends on which is the lightest heavy 
neutrino and the whole transformation E,rVr is needed. Therefore it is convenient 
to work in the physical basis of the heavy neutrinos and to express all the quantities 
in this basis as done in section [7| 
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